gauss forward interpolation formula